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ASYMPTOTIC  GROWTH  OF  CUMULATIVE  AND  REGENERATIVE 
REAM  BREAK-UP  INSTABILITIES  IN  ACCELERATORS 

The  beam  break-up  instability  (BBU)  continues  to  be  a  critical  factor 
which  places  a  limit  on  the  current  and  on  the  pulse  length  in  both  rf  and 
induction  accelerators-  ~  A  transverse  displacement  of  the  particle  beam 
excites  a  non-axisymmetric  mode  in  the  accelerating  structure  (cavity). 

This  non-axisymmetric  mode  further  deflects  the  beam  sideways,  reinforcing 
the  mode  itself.  Depending  on  whether  a  wave  with  negative  group  veloci  ty 
is  present  to  provide  feedback,  BBU  may  either  be  regenerative  or 
cumulative. 

Much  theoretical  effort  on  BBU  in  the  past  twenty  years  has  been 

2  2  s  7  3 

devoted  to  the  cumulative  type,  ’  where  the  accelerating  units  are 

assumed  to  be  decoupled  from  each  other  electromagnetically .  Information 

is  carried  only  by  the  beam.  Under  this  assumption.  Panofsky  and  3ander" 

found  that  the  transverse  displacement  of  the  beam  grows  asymptotically 
1/3 

like  exp(at)  ,  at  a  given  distance  downstream,  when  the  focusing  magnetic 
field  is  absent.  The  model  including  a  general  focusing  magnetic  field  was 
laid  down  in  Ref.  2,  but  it  was  stated  there  that,  except  for  the  case  of 
weak  focusing,  it  is  not  possible  to  obtain  the  asymptotic  growth 
analytically.  In  a  pioneering  paper  published  somewhat  later,  Neil,  Hail 
and  Cooper"’  used  an  entirely  different  approach  and  found  that  the 

1/2 

asymptotic  growth  of  the  cumulative  BBU  behaves  instead  like  exp(bt) 

in  the  presence  of  a  strong  solenoidal  magnetic  field.  Here,  a,  b  are 

parameters  proportional  to  the  beam  current.  These  peculiar  time 

dependences,  at  first  sight,  are  not  expected  from  the  usual  experiences 

with  beam-circuit  interaction.  However,  the  asymptotic  growth  is  firmly 

established,  at  least  in  the  case  of  zero  focusing  magnetic  field,  both  in 

the  "continuum"  limi'-,  (where  coupled  partial  differential  equations  are 
o  3 

solved"’  ),  and  in  the  discrete  model  (vnere  the  individual  beam-cavitv 
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interaction  is  passed  onto  the  subsequent  ones  through  multiplication  of 
5  8 

matrices  ’  ).  From  the  above  asymptotic  growths,  scaling  lavs  on  the  beam 
current  were  established. 

In  this  paper,  we  use  the  continuum  model  and  adopt  a  mode  coupling 

analysis.  The  asymptotic  growth  is  calculated  analytically,  for  both 

cumulative  and  regenerative  3BU,  in  the  presence  of  a  general  focusing 

magnetic  field.  This  work  was  motivated  by  an  attempt  to  understand  tne 

origin  of  the  asymptotic  dependences  mentioned  above,  and  by  the  need  to 

9 

assess  the  importance  of  the  BBU  in  the  two  beam  accelerator  concept 

currently  explored  at  the  Naval  Research  Laboratory.  The  present  analytic 

theory  yields  three  specific  results  which  hitherto  were  not  given  in  the 

l/o 

literature.  First,  the  asymptotic  growth  exp(bt)  exhibited  in  the 

cumulative  BBU  in  the  presence  of  a  strong  solenoidal  magnetic  field^  is  a 

result  of  the  coupling  between  the  slow  beam-cyclotron  mode  and  the  cavity 

1/3 

mode.  Second,  this  growth  is  reduced  to  exp(at)  as  t  ■*  ®.  That  is,  the 

asymptotic  growth  of  the  cumulative  BBU  is  independent  of  the  focusing 

magnetic  field  -  as  if  the  focusing  magnetic  field  were  absent.  Third,  the 

treatment  of  the  cumulative  BBU  is  extended  to  the  regenerative  type  with 

the  inclusion  of  a  negative  group  velocity  v  .  The  exponentiation  factor 

S 

is  modified  by  a  quantity  which  depends  only  on  the  group  velocity,  but  is 
independent  of  the  other  properties  of  the  structure.  Here,  we  shall 
present  the  model  and  the  results.  The  implication  will  be  discussed  but 
the  details  will  be  given  elsewhere. 

For  simplicity,  consider  a  continuous  beam  with  coasting  velocity  v, 
relativistic  mass  factor  y,  and  current  I  streaming  in  a  focusing  magnetic 
field  B  inside  a  series  of  identical  accelerating  units.  Let  c.  (z.t)  be 
the  transverse  displacement  of  the  beam  from  the  axis.  ]  it. r)  be  a  measure 
of  the  deflecting  force  produced  by  the  non-axisymmetric  mode  (with  e1U)ot 
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i. 


I 
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In  the  continuum 


dependence)  in  the  individual  accelerating  units, 
treatment  of  the  cumulative  BBU,  the  governing  equations  for  £  and  q  may  be 
written  as2,3,10 


8,-k)  [»$*’£)] 


+  Y«  S  =  q( z >  t ) 


(1) 


3q(z,t) 

3t 


io)Qq(z,  t) 


“irwo 
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sl(z,  t), 


(2) 


where  a>c  =  je|B/ymoc  and  £  is  the  dimensionless  coupling  constant 
proportional  to  the  beam  current.^  As  usual,  Eq.  (1)  expresses  the 
deflection  of  the  beam  by  the  mode,  whereas  Eq.  (2)  describes  the 
excitation  of  the  mode  by  the  beam's  transverse  displacement  l. 

Assuming  a  dependence  exp(iwt  -  ikz)  for  the  solutions,  Eqs.  (1)  and 
(2)  yield  the  dispersion  relation 


D(«,k)  a 


(3) 


12  13 

In  the  terminology  of  mode  coupling,  ’  this  dispersion  relation 

describes  the  interaction  between  the  cavity  mode  (go  =  aiQ)  and  the  fast  and 

slow  beam-cyclotron  mode  (co  -  kv  =  too^). 

From  the  general  stability  theory,  the  dominant  asymptotic  growth  of 

12 

disturbances  may  be  determined  from  the  Green's  function 


G(z,  t)  =  jda>jdke1“t_lkz  /D(k,«)  -  |d«e1“t_lk(w)z ,  (4) 

r  --  r 

where  the  Bromvi tch  contour  T  lies  sufficiently  far  in  the  lower  half 
complex  oo  plane,  and  k(oo)  is  the  meaningful  solution  obtained  from  the 
dispersion  relation  D(co,k)  =  0. 


Let  us  first  recover  the  previously  established  results  from  Eq.  (4). 

When  the  focusing  magnetic  field  is  absent,  =  0  and 

1/2 

k(co)  a  m/v  +  (coq/v)  [ <*»o e/ ( o>—  <*>q )  J  .  Substitution  of  this  k(w)  in  (4) 

2 

yields  the  following  asymptotic  formula  (from  a  saddle  point  calculation  ): 
|Gn(z,t)  |  -  exp  (l. 64V  1/3),  (5) 

where 

W  a  e  («o3z2/v2}  (t  -  z/v),  (t  >  z/v  >  0).  (6) 

The  asymptotic  solution  (5)  was  first  obtained  by  Panofsky  and  Bander  [cf. 

3 

Eq.  (33)  of  Ref.  (2)]  and  was  confirmed  by  Neil  and  Cooper  and  by 

Gluckstern  et  al.2  in  two  rather  different  analyses.  In  the  other  limit, 

where  a  strong  focusing  magnetic  field  is  present,  the  dominant  interaction 

is  expected  to  be  between  the  (positive  energy)  cavity  mode 

13  14 

(»  -  goo)  and  the  negative  energy  beam-cyclotron  mode,  ’  for  which 
co  -  kv  *  -<oc>  In  this  case,  the  dispersion  relation  (3)  may  be 
approximated  by  (co  -  kv  +  <oc)(co-coq)  =  -coq  J  e/2u>c ,  yielding 

k(co)  =  (w-coc)/v  f  scoq"7(2vcoc(c»>-coo)  ) .  Substituting  this  k(co)  in  Eq.  (4)  and 
performing  a  saddle  point  calculation  similar  to  that  given  in  Ref.  2,  one 
obtains  the  asymptotic  solution 

|G(  z,  t )  |  =  |Gs(z,t)|  -  exp[(2PW/z)1/2],  (7) 

where  p  =  v/co__  and  W  is  given  by  Eq .  (6).  With  the  use  of  the  appropriate 
coupling  constant^  in  W,  the  asymptotic  formula  (7)  is  easily  shown  to  be 
identical  to  the  growth  factor  given  in  Eq.  (5.13)  of  Neil  et  al.3  The 

simple  interpretation,  in  terms  of  mode  coupling,  of  the  asymptotic 

i  /n 

dependence  exp(bt)  evident  in  (7)  is  given  here  tor  the  first  rime.  It 
is  not  readily  extracted  from  the  original  analyses  in  Ref.  5. 


For  general  values  of  to  ,  Eq.  (3)  gives 


1  o  C.W 

,  „  v  a)  1  2  o 

k(0>)  =  -  +  -  tt>  +  - r  , 

V  '  V  -  V  C  ((*>-0)  ) 

0  . 


and  the  saddle  point  contribution  may  also  be  calculated  analytically 
(contrary  to  the  assertion  of  Ref.  2,  at  least  for  a  coasting  beam).  The 
dominant  contribution  to  (4)  gives 


|G(z,  t)  |  -  exp  <J  Re 


1ZW  f  . 

c  _  ,  1 

-  &)  X  1  +  - T— X 

vs  .  3  2 

(  2o>  T 

s 


where  the  dimensionless  time  t  is 


T  =  V(p/z)J,  (10) 

and  u>s  is  the  root  of  the  fourth  degree  polynomial: 

“s3  (1+Ws)  =  <1/2'r)2*  (11) 

It  is  easy  to  show  that  there  is  one  and  only  one  root  of  in  Sq.  'll) 
with  Im  a>s  <  0  for  all  values  of  x  t  0,  and  we  should  use  that  root  of  u>_ 
in  (9). 

The  solution  (9)  implies  that,  given  a  focusing  magnetic  field,  the 
asymptotic  growth  is  independent  of  the  strength  of  the  magnetic  field.  To 
see  this,  consider  a  time  long  enough  so  that  2x  >>  1.  Then  Eq.  (11)  gives 
a>s  a  (1/2t)  e  and  the  solution  (9)  reduces  to  (5),  the  formula 

corresponding  to  zero  focusing  magnetic  field.  This  is  a  rather  surprising 

n 

result,  obtained  directly  from  the  model  of  Panofsky  and  Bander. “  but  is. 
at  first  sight,  contradictory  to  the  findings  of  Neil  et  al." 

The  above  paradox  may  be  resolved  by  noting  that  in  ‘he  "strong 

focusing"  regime,  defined  as  2t  <<  1,  Eq.  (il)  gives 

=  -  i(l/2t)^/,2  [  1- i  ( 1/4)  (2t)  "l/ ^  ] .  Substitution  of  this  .*>  into  (9) 
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yields  (7)  in  the  limit  2x  <<  1.  In  fact,  the  condition  2x  <<  1  may  easily 
be  shown  to  be  similar  to  the  one  imposed  by  Neil  et  alJ  when  they  derived 
Eq.  (7)  using  an  entirely  different  approach.  Thus,  the  dimensionless  time 
X  -  1/2  marks  the  transition  from  the  strong  focusing  to  the  weak  focusing 
regime  [cf  Eq.  (11)].  From  Eqs.  (10)  and  (6),  the  time  (t^)  about  which 

"5 

this  transition  occurs  is  t^,  =  (z/v)(ooc/(Oo)'J/2s. 

As  an  example,  consider  a  1  kA,  30  Me7  electron  beam  coasting  in  a 

solenoidal  magnetic  field  of  3  kG  inside  a  disk-loaded  waveguide  whose 

individual  cavities  support  a  deflecting  TM1  mode  with  oo^/ 2 rt  =  1  GHz. 

Suppose  that  a  misalignment,  say,  produces  an  "excitation"  at  some  axial 

location  and  we  examine  the  BBU  in  response  to  this  excitation  3m 

-4 

downstream.  For  these  parameters,  £  =  4.13  x  10  [cf.  Ref.  11]  and 
tT=  33  ns.  Thus,  if  the  pulse  length  substantially  exceeds  t^,  the 
cumulative  BBU  would  evolve  at  the  later  stage  as  if  the  focusing  magnetic 

field  is  absent. ^ 

The  present  model  may  readily  be  generalized  to  treat  the  regenerative 
3BU  with  the  inclusion  of  a  negative  group  velocity  v  .  The  asymptotic 
formulas  given  below  exolici tly  show  the  change  of  the  character  in  the  BBU 
growth  with  the  sign  of  v  .  The  simplest  way  to  include  the  effects  of  a 
non-zero  group  velocity  is  to  replace  the  factor  in  Eq.  (3)  by 

( )  -  v^(k-kQ),  where  ('joQ,ko)  may  now  be  taken  as  the  point  of 
intersection  of  the  dispersion  curves  in  the  (w.k)  plane  between  the  "beam 
line"  and  that  of  the  slow  wave  structure  formed  bv  the  accelerating  units. 
Vith  this  replacement,  the  Green's  function  (4)  may  again  be  re-evaluated. 
In  the  case  of  no  focusing  magnetic  field,  it  gives 


whereas  in  the  case  of  strong  focusing,  it  reads 


3s(Z’0|  -  exp  |(2pV/z)1/2  (l-vgt/z) 

'  8 


1/2' 


(13) 


where  0  a  v  /v.  It  is  obvious  that  (12)  and  (13)  reduce  to  (5)  and  (7), 

8  8 

respectively,  as  v  ■*  0.  Since  v  <0  for  regenerative  BBU,  it  is  easily 
3  3 

seen  from  Eqs.  (12)  and  (13)  that  both  Gn  and  Gs  grow  like  a  simple 

exponential  function  of  time  [i.e.,  exp(Ct)]  as  t  •*  ®  when  v^  <  0.  This, 

of  course,  is  consistent  with  what  is  expected  from  the  outset  when  a 

12  13 

backward  wave  interacts  with  an  electron  beam  ’  .  It  also  reaffirms  the 

potential  danger  of  the  regenerative  BBU,  as  exp(Ct)  grows  considerably 

1/3  1/2 

faster  than  either  exp(at)  or  axp(bt)  for  large  t. 

The  following  point  may  also  be  of  interest.  A  comparison  of  (12) 

with  (5)  and  (13)  with  (7)  suggests  that  the  modification  in  the 

exponentiation  due  to  a  non-zero  group  velocity  depends  only  on  v  and  is 

5 

otherwise  independent  of  the  accelerating  structure.  It  is  expected,  then, 
that  the  same  modification  would  emerge  in  a  matrix  formulation  similar  to 
the  ones  given  in  Refs.  5  and  8.  The  asymptotic  formulas  (12)  and  (13) 
would  provide  an  immediate  determination  of  the  potential  prevalence  ot  the 
accumulative  or  regenerative  BBU  from  a  knowledge  of  the  structure  (e.g.. 
Brillouin  diagram^’  ^) . 

Finally,  we  remark  that  Eqs.  (12)  and  (13)  are  also  valid  when  0  >  0. 

8 

12 

In  that  case,  the  BBU  becomes  convective  and  there  would  be  an  additional 

restriction  on  the  solutions  (12)  and  (13);  namely,  disturbance  growth  is 

possible  only  for  z/v  <  t  <  z/v  at  a  given  position  z.  In  fact,  all  ot 

8 

the  various  asymptotic  time  dependences  examined  in  this  paper  may  easily 

be  understood  in  terms  of  the  anaiyticitv  ot  the  Laplace  transtorm  ot  t  he 

12 

Green's  function.  From  Briggs'  stability  theory,  it  may  be  shown  that 

the  0  >  0  and  the  0  <  0  cases  correspond  to,  respectively,  anaiyticitv  in 

8  3 


the  lower  half  u>  plane  including  the  real  u  axis,  and  a  branch  point 


somewhere  in  the  lower  half  <o  plane.  The  special  case  0  =0  (cumulative 

s 

BBU)  corresponds  to  analyticity  in  the  lower  w  plane,  but  with  a  branch 
point  on  the  real  o>  axis,  leading  to  an  exponential  growth  whose  exponent 
is  a  fractional  power  of  time,  as  exhibited  in  Eqs.  (5)  and  (7). 

In  addition  to  the  general  theory  reported  above,  a  similar 
formulation  has  been  carried  out  to  calculate  the  excitation  of  the  BBU  in 

the  accelerated  beam,  by  the  misalignments  in  the  primary  beam,  in  the  two 

q 

beam  accelerator  configuration  recently  proposed  by  Friedman  and  Serlin. 
The  preliminary  conclusion  is  that,  in  the  parameter  regime  currently 
studied,  the  BBU  does  not  appear  to  be  serious. 

I  would  like  to  thank  M.  Friedman  and  D.  Chernin  for  helpful 
discussion.  This  work  is  supported  by  the  U.  S.  Department  of  Energy, 
Contract  No.  QE-AI05-36-ER13585 . 
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different  regimes,  as  the  stabilizing  influence  due  to  a  finite  0  may 
be  accounted  for  phenomenologically  by  simply  multiplying  the  Green's 
function  by  exp  ( -<a^~/2Q) . 

11.  The  coupling  constant  z  depends  on  the  accelerating  structure  and 
on  the  deflecting  mode  under  consideration.  The  configuration 
treated  in  Ref.  3  consists  of  a  circular  waveguide  loaded  with 
identical  apertured  disks  aiong  rhe  guide  axis.  In  that  case, 

z  =  0.0248  (I/lkA)6/r  if  the  deflecting  mode  within  the  individual 
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In  the  more  general  configuration  treated  in  Refs.  2  and  5.  the 
accelerating  units  are  separated  by  a  distance  L  and  the  n-th  unit  is 

->  2 

located  at  z  =  nL.  In  that  case,  e  =  (v  K/a)Q  L)  I/(170ykA)  where  k 
is  proportional  to  the  "transverse  impedance"  of  the  structure. 

[ K  has  a  unit  of  inverse  length;  it  is  identical  to  the  It  defined  in 
Eq.  (3.11)  of  Ref.  5j.  Note  that  £  is  independent  of  the  focusing 
magnetic  field  but  is  inversely  proportional  to  y. 

R.  J.  Briggs,  Electron  Stream  Interaction  in  Plasmas,  MIT  Press, 
Cambridge,  MA,  (1964),  Ch.  2. 

See  e.g.,  M.  V.  Chodorow  and  C.  Susskind,  Fundamentals  of  Microwave 
Electronics .  (McGraw-Hill,  New  York,  1964)  tor  a  general 
discussion  of  mode  coupling  and  beam-circuit  interaction. 

A.  Bers  and  S.  Gruber,  Appl.  Phys.  Lett.  6,  27  (1965);  L.  S.  Hall  and 
'7.  Heckrotte.  Phys.  Fluids  9-  ’-496  (1965). 

It  should  be  emphasized  that  the  damping  mechanisms  associated  with 
phase  mixing  and  a  finite  Q  might  have  a  decisive  influence  on 
the  3BU  growth.  Nevertheless,  from  an  analytical  point  of  view, 
it  is  interesting  that  the  asymptotic  growth  in  an  unstable 
system  described  by  the  dispersion  relation  (3)  is  independent  of 
In  the  :ase  of  an  accelerating  beam,  a  careful  examination  of  the 
Green's  function  would  show  that  the  asymptotic  growth  of  the 
cumulative  BBU  is  also  independent  of  the  focusing  magnetic  field, 
in  the  sense  described  here.  In  fact,  the  analyses  and  rhe 
conclusions  given  in  this  paper  need  not  be  restricted  to  a  solenoidal 
focusing  magnetic  field.  For  general  linear  transverse  focusing,  we 
simDly  replace  bv  k^c,  where  k^  is  the  betatron  wave  number,  ind 

the  asvmDtotic  formulas  remain  valid. 
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